Evidence for strong electronic correlations in the spectra of Sr 2 Ru04 
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The importance of electronic correlation effects in the layered perovskite Sr2Ru04 is evidenced. 
To this end we use state-of-the-art LDA+DMFT (Local Density Approximation + Dynamical Mean- 
Field Theory) in the basis of Wannier functions to compute spectral functions and the quasiparticle 
dispersion of Sr2Ru04. The spectra are found to be in good agreement with various spectroscopic 
experiments. We also calculate the k-dependence of the quasiparticle bands and compare the re- 
sults with new angle resolved photoemission (ARPES) data. Two typical manifestations of strong 
Coulomb correlations are revealed: (i) the calculated quasiparticle mass enhancement of m*/m ~ 2.5 
agrees with various experimental results, and (ii) the satellite structure at about 3 eV binding energy 
observed in photoemission experiments is shown to be the lower Hubbard band. For these reasons 
Sr2Ru04 is identified as a strongly correlated Ad electron material. 

PACS numbers: 71.27. +a, 79.60.-i 

I. INTRODUCTION 

Intensive research on Sr2Ru04 began after the discovery of superconductivity at temperatures below lKi Since it 
is widely believed that this system may help to clarify the mechanism behind high-T c superconductivity, considerable 
theoretical and experimental effort was put into the investigation of this material which has unconventional properties 
in spite of a relatively simple electronic structure. While it is generally recognized that in 3d transition metal 
compounds electron correlations play a crucial role^ the question whether the Ad system Sr 2 Ru04 should be considered 
a strongly correlated system, too, or whether its electronic properties can be understood within conventional band 
theory remained open. Indeed, the 4c? states of the Ru-ion are more extended than the 3d states and hence correlation 
effects should be less significant than, e.g., in high-T c cuprates. On the other hand, the effective quasiparticle mass 
obtained from de Haas- van Alphen (dHvA)^ ARPES^ and infrared optical experimental is 3-4 times larger than 
the results obtained from standard band calculations*!*^ The temperature-independent contribution to the magnetic 
susceptibility and specific heat constant are also significantly larger than that given by LDA£ These facts indicate 
that electron correlations play an important role in Sr2Ru04. 

During the last decade intensive studies of Sr2Ru04 using various spectroscopic techniques were 
pcrformediSiifliiLiSiiiiii^ The evidence derived from these experiments which point towards correlation effects can 
be summarized as follows: (i) the measured bandwidth and density of states (DOS) at the Fermi level deviate by 
roughly a factor of 32ii£ from bandstructure calculations^**^ and (ii) there is a peculiar satellite at —3 eV in the 
photoemission spectra (PES)i 12 i 13 i 14 i 15 

The presence of a satellite structure in PES on transition metal oxides, now commonly interpreted as the lower 
Hubbard band (LHB), is generally taken as strong evidence for the importance of correlations. Such a satellite was 
first experimentally observed by Fujimori et al*^ for the d 1 perovskite-type Ti 3+ and V 4+ oxides. By applying the ah 
initio LDA+DMFT approachpiiSiiSiSSiSiiSS the corresponding structure in the many-body spectrum was later indeed 
identified as the LHBi 1 ^* 1 ^^*^* 2 ^ While in d 1 compounds the 3d-band is usually well separated from the oxygen 2p 
band (making the experimental observation and theoretical interpretation as the LHB relatively simple) the energy 
separation between the Ru-4d and 0-2p states in Sr2Ru04 is much smaller (see below). In this case the LHB may 
overlap with the oxygen 2p bands, making the interpretation of structures in the experimental spectra ambiguous. The 
difference in experimental conditions (photon energy, surface sensitivity, sample and surface quality, etc) complicate 
the situation even more. Therefore interpretations of spectroscopic data for Sr 2 Ru04 are often controversial. 

Early investigations of the electronic structure of Sr2Ru04 used the local density approximation (LDA) to reveal 
the similarities and differences with the electronic properties of cupratc superconductors*!*^ It was then proposed 
that the superconductivity of Sr 2 Ru04 may be unconventional, namely of triplet type(*2£*2£ see the comprehensive 
reviews *S**2*i A quantitative model for triplet superconductivity based on first principles calculations for the electronic 
structure and magnetic susceptibility was suggested by Mazin and Singh*S The electronic structure of Sr2Ru04 and 
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Sr2Rh04 was compared in* and the possibility of a magnetic ground state of Si'2Ru04 was studied within the general 
gradient approximation (GGA) — 

The Fermi surface of Sr2Ru04 was also investigated by LDA. According to these studies, the Fermi surface consists 
of three cylindrical sheetsjSiS* 3 ^ in agreement with dHvA experiments^ By contrast, ARPES experiments predicted 
a significantly different Fermi-surface topology— In principle, such a discrepancy may be due to strong elec- 
tronic correlations which are not taken into account in LDA. However, detailed photoemission studies^ and scanning- 
tunneling microscopy^ subsequently discovered a surface reconstruction which seemed to resolve the controversy^Zi3& 

The importance of correlation effects was studied by Perez-Navarro et a/— and Arita et a/.— Although correlations 
are not very strong, their inclusion was found to be important for a proper description of the electronic structure 3 ^ 
and a microscopic understanding of superconductivity— 

A first multi-band investigation based on the dynamical mean-field theory (DMFT) to clarify the discrepancy 
between dHvA 3 . and photoemission 3 ^ 3 ^ data was performed by Liebsch and Lichtenstein— They observed a charge 
flow from the narrow xz, yz bands to the wide xy band leading to a shift of the van Hove singularity close to Ep, 
and derived quasiparticle bands by self-consistent second-order perturbation theory for the self-energy, finding a mass 
renormalization of 2.1-2.6— in agreement with experiments^ Anisimov et alm& investigated the isoelectric series of 
alloys Ca 2 _ x Sr x Ru0 4 by means of LDA+U for x = and DMFT(NCA) for 0.5 < x < 2.0. In the latter doping range 
the scenario of an orbital selective Mott transition (OSMT) was proposed. 

In this paper we address the question of the importance of correlation effects in Si"2Ru04 by a realistic 
LDA+DMFT(QMC) calculation within a Wannier function (WF) formalismi 3 . This improved LDA+DMFT(QMC) 
scheme allows one to take into account the influence of correlated orbitals (4d-t2 fl orbitals of Ru in our case) on 
all other states. This is essential when d and oxygen p states overlap as is the case in Sr 2 Ru04. The comparison 
between our theoretical results and experiment clearly shows that electronic correlations have a strong influence on 
the electronic structure and lead to the formation of a pronounced LHB in the spectrum. 

The paper is organized as follows. In scction[H]we present results for the band structure obtained by LDA (subsection 
II. A) and LDA+DMFT(QMC) in the WF basis (subsection II. B), respectively. Section UTTI contains a comparison of 
our LDA+DMFT(QMC) results with XPS (x-ray photoemission spectroscopy), XAS (x-ray absorption spectroscopy), 
and NEXAFS (near edge x-ray fine structure) experiments, as well as with new experimental PES (photoemission 
experiments) data IjHI A|) and recent ARPES experiments (subsections IIII Bl and IIIl"Cll . We conclude the paper with 
a summary. I1VI 

II. ELECTRONIC STRUCTURE 
A. LDA band structure 

Sr2Ru0 4 has the undistorted single-layered K2NiF4-type structure with the space group IA/mmm and lattice 
parameters a=&=3.8603 A, c=12.729 A— The structure is formed by layers of RuOg-octahedra separated by Sr-ions. 
The RuOg-octahedra are slightly elongated along the c-axis. Therefore the coordination of Ru-ions locally has a 
tetragonal symmetry. 

Our first- principle calculation of the electronic structure of Sr2Ru04 is based on density functional theory (DFT) 
within the LDA approximation's^ using the linearized muffin-tin orbitals (LMTO) method— The partial densities 
of states for Sr 2 Ru04 are shown in Fig. ^ They are in good agreement with results of previous calculationsi 6 i 7 i 8 i 31 i 32 
The strontium Ad states are almost empty and lie above 3 eV; the 0-2p derived bands are filled and extend from -8 
eV to -1 cV. 

Physically most interesting are the partially filled ruthenium 4<i-statcs. Due to the octahedral coordination of the 
oxygen ions, the Ru-4d states are split into t2 g and e g orbitals (see Fig.[2J). Owing to the stronger hybridization of 
the two e g orbitals with oxygen p-states the corresponding bands lie above the three t 2g bands in the energy region 
from 0.5 eV to 5 eV. In Sr 2 Ru04 four Ru-4d electrons occupy the three t 2g bands (d 4 configuration). The layered 
crystal structure of Sr2Ru04 results in a two-dimensional DOS of the xy orbital while the xz, yz orbitals have nearly 
one-dimensional character (sec Fig.|2}. The xy orbital hybridizes with xy orbitals of the four Ru neighbors and thus 
has a bandwidth almost twice as large (W xy =2.8 eV) as that of the xz, yz orbitals (W XZi y 2 =1.5 eV) which hybridize 
with corresponding orbitals of two Ru neighbors only. 

The LDA bands in the energy window —3 ... 1 eV are shown in Fig. |3J In contrast to a typical d 1 systems^ there 
is no well-pronounced separation of the oxygen 2p and ruthenium Ad states in Sr2Ru04. More precisely, Fig. [3] shows 
that while the Ru-4<i xz, yz orbitals are separated from the oxygen 2p bands the Ru-4<i xy orbital strongly overlaps 
with these oxygen bands. 
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FIG. 1: Partial LDA DOS for SnRuCU. The Fermi level corresponds to zero. 
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FIG. 2: Orbitally projected LDA Ru-4d DOS. The Fermi level corresponds to zero. 



B. LDA+DMFT(QMC) results: effect of correlations 

Sr2Ru04 is a paramagnetic metali It is well known that the paramagnetic state of a correlated metal is well 
described by the DMFT (for reviews see^^Si). Within DMFT the electronic lattice problem is mapped onto a 
single-impurity Anderson model with a self-consistency conditional^ This mapping, which becomes exact in the 
limit of large coordination number of the lattice^ allows one to investigate the dynamics of correlated lattice elec- 
trons non-perturbatively at all interaction strengths. We use the LDA+DMFT ab initio techniqueiiiSii2i2^ (for an 



4 




FIG. 3: (Colour online) S^RuCU LDA band structure along high symmetry directions of the Brillouin zone. Light line - LDA 
bands, dark line - bands obtained by Wannier function projection on t2 g orbitals. The Fermi level corresponds to zero. 

introduction see^Si for reviews se e22i2L££i££) to investigate correlation effects in Sr2Ru04. The effective impurity prob- 
lem corresponding to the many-body Hamiltonian is solved by quantum Monte Carlo simulations^ The LDA+DMFT 
approach was recently improved by employing a Wannier functions (WF) formalism^ which allows one to project 
the Hamilton matrix from the full-orbital space to a selected set of relevant orbitals. The projection ensures that the 
information about all states in the system is kept. In the present work we use the WF formalism to construct an 
effective few-orbital Hamiltonian with t2 g symmetry and to take into account the influence of correlated t2 9 -orbitals on 
other states. A three-orbital Hamiltonian obtained by the WF projection with dispersions presented in Fig. [3] (black 
lines) was used as an ab initio setup of the correlation problem. Ab initio values of the orbitally averaged Coulomb 
interaction parameter [7=1.7 eV and Hund exchange parameter J=0.7 eV were obtained by constrained LDA cal- 
culationsS We emphasize that not only the on-site e g screening and the screening from Ru-ions of the RuC>2 plane 
were taken into account in the calculation of the Coulomb interaction parameter but also screening from neighboring 
Ru02 planes. 

In the particular case of three t2 9 -orbitals U is equal to the inter-orbital Coulomb repulsion U'M^ Thus we 
obtain U = U' + 2 J =3.1 eV for the local intra-orbital Coulomb repulsion. We note that our values of U and J differ 
substantially from those by Liebsch and Lichtenstein^i who assumed a much smaller Hund exchange parameter. These 
authors estimated the Coulomb repulsion parameter from the XPS spectrum^ using the position of the resonance 
satellite. The value of about 1.5 eV obtained thereby agrees well with our calculated value for U. 

The three-orbital, projected Hamiltonian together with the ab initio Coulomb interaction parameters were used as 
input for the QMC simulation of the effective quantum impurity problem arising in the DMFT. The simulations were 
performed for an inverse temperature /3=10 eV -1 using 40 imaginary time slices (Ar=0.25). The imaginary time 
QMC data were analytically continued by maximum entropy^ The results are shown in Fig. We find a pronounced 
lower Hubbard band (LHB) between -5 and -1 eV, a quasiparticlc peak (QP) around the Fermi level, and an upper 
Hubbard band (UHB) at about 1.5 eV. Real and imaginary parts of the corresponding self-energy for real frequencies 
(for details see Appendix B in££) for t2 9 orbitals are shown in Fig. |SJ The mass enhancement calculated from the 
derivative of ReE at the Fermi level amounts to 2.62 for the xy orbital and 2.28 for the xz, yz orbitals, in agreement 
with results from ARPES^ dHvA^ and infrared optical experiments^ A detailed analysis of the structures seen in 
Fig. 2] is presented in Appendix B 

After having calculated the self-energy for real frequencies one may perform the inverse transformation from 
the reduced Wannier basis back to the full LMTO basis^ii This step allows one to take into account the influence 
of the three correlated t2 S orbitals on all other states of ruthenium, oxygen and strontium. The comparison of the 
noninteracting LDA partial density of states with the one obtained by the inverse transformation is shown in Fig. [5] 
Since the hybridization of Ru and O is quite strong the oxygen states are changed rather significantly by correlation 
effects on Ru-ions. These changes are most pronounced in the energy region between -4 and -1 eV. One can see 
that apical oxygen atoms are more affected by correlations than in-plane atoms. We believe that this is due to the 
one-dimensional character of the xz,yz states of the Ru-4d shell. Hence correlation effects are much stronger for 
these orbitals; consequently the DOS for the apical oxygen atoms is strongly modified. The strontium states are less 
affected. 

A comparison between the partial LDA DOS of Ru-4g? and the DOS obtained using the full-orbital self-energy from 
our LDA+DMFT(QMC) calculations is shown in Fig.0 The main effect of the correlations on the Ru site is seen to 
be a transfer of spectral weight from the energy region near the Fermi level to the lower and upper Hubbard bands 
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FIG. 4: (Colour online) Ru-4d(t,2 g ) spectral functions obtained within LDA+DMFT (QMC) using a projected Hamiltonian. 
Dark curve: xy orbital, light curve: xz, yz orbitals; £7=1.7 eV, J=0.7 eV. The Fermi level corresponds to zero. 




FIG. 5: (Colour online) Self-energy on the real energy axis for xy (black line) and xz, yz (light line) orbitals of Sr2Ru04. Solid 
line: real part; dashed line: imaginary part. The Fermi level corresponds to zero. 



range from -4 eV to -1 cV, and from 1 eV to 2 eV respectively. 



III. COMPARISON WITH EXPERIMENTAL DATA 



A. XPS and PES experiments 



We will now compare the computed LDA+DMFT(QMC) spectral functions for the t2 g electrons (light solid lines in 
Fig. ISll§l llL)fl and those calculated using the full-orbital self-energy (black solid lines in the same figures) with several 
experimentally obtained spectra describing both valence and conduction bands. To compare with experiment we took 
into account the photoemission cross section ratio for Ru-4<i and 0-2p states as a function of photon energyiS We 
found that, in general, an energy dependent broadening of the theoretical spectral functions gives better agreement 
with the experimental data (see££ and^i). For this the theory curves were convoluted using a Gaussian with a full 
width at half maximum increasing as C-E+g. Here E is the binding energy, g is the experimental resolution, and 
C characterizes the increase of the broadening with energy upon moving away from Fermi level due to corc-holc life 
time effects. The maximally allowed broadening was restricted to 1 cV. Specific values of C and g parameters used 
for comparison with experiment are indicated in the corresponding figures. 
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FIG. 6: (Colour online) Comparison of the total and partial LDA DOS (light curve) and the DOS calculated using the 
full-orbital self-energy from LDA+DMFT(QMC) (black curve). The Fermi level corresponds to zero. 



0.5- 



I 0.5 
o 

> 

S3 

Bo.5 
o 

Q 





0.5 




Ru-4J 


t 2g (xy) 

\-^L .1,1. 




t 2g fe yz) 
.1,1. 




2 2 

e g (3z -r ) 


| LDA 

-1 — LDA+DMFT(QMC) ( 

1 T*^ - 1 1 . 


e g (x -y ) 

full-orbital) 



-8 -6 -4 -2 2 4 6 8 
Energy (eV) 



FIG. 7: (Colour online) The same as in Fig. but for Ru-4d states only. 
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1. Comparison with previous XPS experiments 
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FIG. 8: (Colour online) Theoretical spectral functions of Sr 2 Ru0 4 calculated by LDA+DMFT(QMC); light solid line: t 2g 
orbitals, black solid line: full-orbital self-energy. In (a) the results are compared to an angle-integrated valence band photoe- 
mission spectrum— obtained with a photon energy E p h = 30 eV, while in (b) we compare to an XPS spectrum obtained with 
a photon energy E p ^ = 1486.6 eV^- The theoretical spectra are convoluted using a linear broadening -0.04-E+0.25 for (a), 
and -0.14-E+0.25 for (b) to account for the experimental resolution. Intensities are normalized on the area under curves. The 
Fermi level corresponds to zero. 



In Fig. [S] we compare an angle integrated valence band photo-emission spectrum of S^RuO^ii [Fig. Eta)] and the 
XPS spectra from Refii^ [Fig. Efb)] with the theoretical spectral functions. The contributions from Ru-4d and 0-2p 
spectra were weighted according to the photocmission cross section ratio^S 3:1 for Fig.|SJa), corresponding to a photon 
energy 30 eV, and 40:1 for Fig. IHtb), corresponding to a photon energy of 1486.6 eV. The theoretical spectra were 
multiplied with the Fermi function at T = 150 K and 300 K, respectively. In Fig.[S]a linear broadening was employed. 
UPS (ultra-violet photoemission) data of the valence band of Si^RuO^ obtained at higher photon energy of 60 eV 
and 110 eV show similar features as the PES and XPS spectra in Fig. [HI 

We note that for energies below -2 eV the weight of oxygen states in the spectrum (according to the cross section 
ratio) becomes essential. Indeed, only by proper inclusion of this contribution can a satisfactory description of the 
experimental spectra with the full-orbital LDA+DMFT(QMC) spectral function be achieved. 

The experimental spectrum in Fig. [HJa) was obtained at a rather low photon energy. Therefore, according to the 
cross section ratio, the contribution of oxygen states is considerable. The calculated curves reproduce all features of the 
complicated structure of the experimental data, although the positions and weight of peaks agree only qualitatively. 
The most serious disagreement with experiment can be observed in the energy region from -2.5 eV up to the Fermi 
level. Despite similar line shapes the theoretical curves have too little weight. This failure can, for example, be 
attributed to an underestimation of the oxygen contribution in the theoretical curves, or to matrix element effects 
which may be significant at low photon energies. Nevertheless, the LDA+DMFT(QMC) spectra are seen to be in 
much better agreement with experiment than the LDA results. In particular, in LDA+DMFT the LDA peak near 
-0.5 eV is merely becomes a plateau - the redistribution of spectral weight obviously being an effect of correlations. 
Comparing the theoretical t2 9 and full-orbital spectra with experiment, one can see that the latter yield a much better 
description in a wide energy range. 

The XPS spectrum in Fig. HJb) obtained at a very high photon energy is seen to be almost exclusively Ru-4d 
states. Obviously, the full-orbital spectral function gives good agreement with the XPS data. Moreover, one observes 
a pronounced maximum at -3 eV which experimentalists previously interpreted as the LHB^S This conjecture is now 
confirmed by our calculations (see the detailed discussion below). 



2. Comparison with new PES experiment 

Clean (001) surfaces of high purity single crystal samples were obtained by in situ cleavage at 20 K in ultrahigh 
vacuum. Angle-integrated and angle-resolved spectra were measured using the GAMMADATA-SCIENTA SES200 
Analyzer at BL25SU of SPring-8 by use of circularly polarized light. Both measurements were performed at 700 eV by 
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FIG. 9: (Colour online) Similar plot as in Fig.[S]but now we compare with PES spectrum of Sr2Ru04. The theoretical spectra 
are convoluted using a linear broadening -0.04-E+0.20 to account for the experimental resolution. 

detecting near normal emission electrons to obtain highest bulk sensitivity at this photon energy. The resolution for 
these measurements was set to 200 meV (Fig. 01 and 120 meV (Fig. lll|) . respectively. The energy scale was calibrated 
by the Fermi edge of Au. The surface cleanliness was checked by the absence of possible additional spectral weight 
on the higher binding energy side of the O Is peak, and by the absence of the C Is contribution. 

In Fig. [5] we compare photoemission spectrum of Sr2Ru04 with spectral functions calculated by 
LDA+DMFT(QMC). A weighted sum of Ru-4d and 0-2p spectral functions according to the photoemission cross 
section ratio 17:1— was used - corresponding to an experimental photon energy of 700 eV. Theoretical spectra were 
multiplied with the Fermi function at 20 K and were linearly broadened to account for the experimental resolution. 



We now discuss and interpret the experimental and theoretical spectra, and also check the presence of a LHB 
in the computed spectra. In Fig. |Hfb) the structure in the experimental XPS-spectrum at -3 eV was interpreted 
as the LHB— In the full-orbital LDA+DMFT(QMC) spectral function (black solid line) a corresponding feature 
is indeed visible, but has less intensity and appears only as a shallow shoulder rather than a distinct bulge. This 
structure is obviously not described by LDA (dashed line). Looking at the light solid line in Fig.[H{b) which represents 
the LDA+DMFT(QMC) DOS for the t2 ff orbitals alone, we are able to identify this shoulder unambiguously as a 
result of the LHB. Thus we have theoretically confirmed the interpretation of Yokoya et al.H A very similar feature 
corresponding to the LHB in our LDA+DMFT(QMC) DOS was reported in Ref^. 

The situation is similar in the PES spectrum. Because of lower photon energies, and due to the enhancement of the 
0-2p contribution in the PES spectrum seen in Fig. EI one again cannot identify the -3 eV satellite directly. However, 
one can recognize a feature in the spectrum whose position coincides with the LHB in our LDA+DMFT(QMC) 
calculations. In the first theoretical DMFT work on the ruthenate^i model-Coulomb parameters were chosen as 
17=0.8 eV and J=0.2 eV. As a consequence a less well-defined LHB with low spectral weight was obtained. 



In Fig. llOf a) the O-ls x-ray absorption spectrum (XAS) of Si^RuO^i^ representing the conduction band is compared 
with the 0-2p spectral function calculated via LDA-fDMFT(QMC). Furthermore, Fig. HOf b) shows the near edge x- 
ray fine structure spectrum (NEXAFS) of S^RuO^ together with the theoretical spectral functions. The theoretical 
spectra are multiplied with the Fermi function at T = 300 K and convoluted using linear broadening to account for 
the experimental resolution. 

The agreement between theory and experiment in Fig. llOl is found to be only qualitative. This may be due to empty 
states (conduction band) in the LMTO method. Namely, the conventional LMTO choice of the MTO linearization 
energy point lies inside the occupied part of the bands. As a result the unoccupied states in LDA calculation are not 



3. Interpretation 



B. XAS and NEXAFS experiments 
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FIG. 10: (Colour online) Theoretical spectra of Sr2Ru04 calculated by LDA+DMFT(QMC); light solid line: only t2 9 electrons, 
black solid line: full-orbital self-energy using 0-2p states, (a) Comparison with O-ls XAS spectrum^- (b) comparison with 
NEXAFS spectrum!- Theoretical spectra are convoluted using linear broadening 0.0TE+0.35 for (a), and 0.02-E+0.2 for (b) 
to account for the experimental resolution. Intensities are normalized on the area under curves. The Fermi level corresponds 
to zero. 



properly described. Nevertheless, due to the spectral weight redistribution in DMFT calculation, the agreement with 
experiment is improved (black solid line) in comparison with the bare LDA 0-2p DOS (light dashed line). 



C. High-energy bulk ARPES experiment 



The comparison of experimental and theoretical spectra presented in the previous sections proves the existence of 
a LHB in Sr2Ru04. Hence Sr2Ru04 must be considered a strongly correlated electron material. Further evidence 
for correlation effects comes from the renormalization of quasiparticlc properties, i.e., effective masses and band 
dispersions. This renormalization was addressed by Liebsch and Lichtenstein-i who derived quasiparticlc bands from 
self-consistent second-order self-energy. The Coulomb repulsion and Hund's rule exchange parameters chosen by these 
authors were rather small ([7=1.2-1.5 eV, J=0.2-0.4 eV), resulting in a LHB with low spectral weight; nevertheless 
the effective masses (m* ~ 2.1 — 2.6) were found to be in good agreement with experiments. In our investigation 
we determined the local Coulomb repulsion by the ab initio constrained LDA yielding [7=1.7 eV and J=0.7 eV, i.e., 
a value of [7=3.1 eV which is twice as large as that used in Ref4i. This value of U now produces a pronounced 
LHB and, at the same time, gives almost the same values of m* /m as those reported in Ref. — : m*/m=2.62 for the 
xy orbital and 2.28 for the xz, yz orbitals, respectively. All these values are in good agreement with experimental 
estimations 

We will now proceed to calculate k-resolved spectra and quasiparticlc band dispersions following the strategy 
proposed in RefM£2£& and employed by us to calculate the ARPES spectra of SrVC^^ ARPES computations have 
previously been performed also for the 2D Hubbard model by Maier et al~ in the framework of the dynamical cluster 
approximation (DCA)^ and by Sadovskii et alm& within the so-called DMFT+Efc approach. 

To compare theoretical quasiparticlc bands with the dispersion extracted from high-energy ARPES data we first 
calculate the fc-resolved spectral function A(k,u) for Sr^RuCU (for details see§£) defined by 

A(k,w) = --ImTrG(fc,w). (1) 

This quantity is determined by the diagonal elements of the Green function matrix in orbital space 

G(fe,a;) = hSH-J/fffe)]- 1 , (2) 

where H WF is the few-orbital Hamiltonian with t2 S symmetry obtained by WF projection. The corresponding 
eigenvalues are pictured in Fig.[3]as black lines. Since QMC only provides the self-energy E for Matsubara frequencies 
iui n and the local spectrum A(uj), the calculation of A(k,u) requires a method to compute E for real frequencies u>. 



10 




FIG. 11: (Colour online) Comparison of the dispersion extracted from high-energy ARPES data& and LDA-fDMFT(QMC). 
Experimental data contain the second derivative of the energy distribution curves, peak positions of the second derivative (black 
dots), and kF estimated from the momentum distribution curves at Ef (triangles). The theoretical dispersion is indicated by 
open circles. 



This is achieved by first employing Kramcrs-Kronig to obtain 

oo 

G(u + irj) = I dJ — . (3) 

J LO - LO 1 + IT) 

— oo 

The local Green function and the complex self-energy are related by the fc-integrated Dyson equation 

G{lu) = fdk [uj + fx- S(w) - hf(k)}- 1 . (4) 



BZ 

Solving Eq. Q for given G(u>) for the self-energy £(w) leads to the results presented in Fig. EI which are then inserted 
into Eq. J5J) to obtain the spectral function A(\t, ui) from Eq. Q. The maxima of A(k, to) obtained from this procedure 
are shown in Fig. llll as black dots without error bars. Compared to experimental results our theoretical quasiparticlc 
bands are shifted by -0.08 eV. 

ARPES data directly provide the energy distribution curvesjSLffi whose second derivative represents the dispersion, 
which is shown in Fig. llll bv a "rainbow scale", red being the highest energy. Closed circles with error bars denote the 
peak positions of the second derivative. The Fermi momenta kp estimated from the momentum distribution curves 
(MDCs)££ at Ep are shown by triangles. We see that in the experiment there are two bands crossing the Fermi level 
at different k-points in the (0,0)-(7r,0) direction. There are also two bands in the calculated dispersion curve but they 
are almost degenerate. Nevertheless the position of the bottom of quasiparticle bands, the intersection with the Fermi 
level, and the shape of the experimental and theoretical dispersions qualitatively agree. 
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IV. CONCLUSION 

In this paper we focused on the long-standing controversy concerning the strength of correlation effects in Si^RuO,!, 
i.e., on the question whether Sr 2 Ru04 should be considered a strongly correlated electron material. This is generally 
the case if the ratio of Coulomb interaction and kinetic energy (bandwidth) is larger than unity. In particular, electronic 
correlations lead to a typical redistribution of spectral weight and thereby to the formation of well pronounced lower 
and upper Hubbard bands (LHB, UHB). We note that even on the experimental level the unambiguous identification 
of a maximum in the spectrum of Sr2Ru04 is made complicated by the overlap of the Ru-4d and 0-2p bands. 

To answer this question we first calculated the electronic structure of Sr 2 Ru04 within the conventional band theory 
using LDA. The correlations were then taken into account in the framework of the LDA+DMFT(QMC) scheme using 
ab initio values for the Coulomb and Hund exchange parameters. We found that the ratio of Coulomb interaction and 
bandwidth is indeed larger than unity in Si - 2Ru04, despite the rather extended 4d-states of Ru, leading to a distinctive 
redistribution of spectral weight and to the formation of a well-pronounced LHB. By comparing our theoretical spectra 
with XPS experiments we unambiguously identified this LHB with the structure observed at -3 eV. By contrast, the 
LDA DOS shows no distinctive feature at that energy. 

To describe the experimental spectra in a wide energy range we employed a Wannier function formalism to transform 
the self-energy operator back to the full-orbital basis. The theoretical spectra obtained in that way agree very well 
with high-photon energy photoemission data; in particular, they reproduce the shoulder in the spectrum caused by 
the LHB. The basic features of the low photon energy UPS and intermediate energy PES spectra are also reproduced 
by the LDA+DMFT(QMC). Quasiparticle bands induced by correlations with mass renormalization of about 2.5 
agree well with results from ARPES, dHvA and infrared optical experiments. The LDA+DMFT(QMC) derived 
quasiparticle bands are even in quantitative agreement with the dispersion extracted from ARPES data. Taken 
together these results provide clear evidence for strong electronic correlation effects in S^RuC^. Hence, although 
Sr 2 Ru04 is a 4c? system it must be regarded as a strongly correlated electron material. 
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APPENDIX: ANALYSIS OF LDA+DMFT(QMC) SPECTRUM 

An important first step in interpreting the structures in the DOS in Fig. 01 is to identify purely atomic excitations. 
These will of course be shifted and broadened by correlation effects, but should nevertheless still be prominent. These 
excitations can readily be obtained from the atomic level picture. We denote the bare level energy of the xz, yz 
orbitals as t X z,yz = eo- Note that we do not know this value a priori, since it is neither the center-of-gravity of the 
bare LDA-DOS nor a particular peak position in the correlated DOS. It is, in fact, the unknown shift due to the 
double-counting correction of the LDA. 

Hund exchange and Coulomb parameters from constrained LDA are J = 0.7 eV and U = 1.7 eV. For three t2 9 - 
orbitals we have V = U = 1.7 eV and thus U = U' + 2 J = 3.1 eV. Finally, the splitting of the centers of gravity of 
the xy and xz, yz DOS is Ae = e xy — e X z,yz = 0.1 eV. 

With these information we can construct the basis spanning the ground state and then calculate the excited states 
and their relative energies. Concentrating on a particular basis vector of the ground state manifold, we obtain the 
scheme in Tab. [I] where we show the states contributing to the possible single-particle excitations. The corresponding 
excitation energies with the unknown level shift eo are listed in the third column. Obviously, the transition c xy \Gs) 
represents the excitation with the lowest energy of the xy-DOS, i.e., should be identified with the position of lowest 
peak at energy E x = -3.2 eV in the LDA+DMFT(QMC) DOS in Fig.H leading to e = ~7 eV and the final numerical 
values in the last column of Tab. [I] 

The identification of the further structures is now straightforward. For the xy orbital, we must attribute the peak 
around E = 0.5 eV to the excited states c\ \GS). Likewise, in the xz, yz manifold the peak at E « —3.5 eV 
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TABLE I: Ground-state as well as single-particle excitations and their energies. 

\xz) \yz) \xy) Energy E a Exitation energy e value, eV 

Gs Tl> t) t) 4e + Ae + 617 - 13 J 

occupied states 
e = 7?gs — S Q 





ITI) 


r> 


|0) 


3e +3(7-57 


e + Ae + 3(7 - 8J=e + 3.8 


-3.2 


C T 


II) 


IT) 


IT) 


3e + Ae + 3(7 - 7J 


e + 3(7 - 6J=e + 5.1 


-1.9 




ITI) 


|o> 


IT) 


3e + Ae + 317 - 5J 


e + 3(7 - 8J=e + 3.7 


-3.3 


c J 


It) 


IT) 


IT) 


3e + Ae + 3U - 9J 
empty states 
e = E a - E G s 


e + 3(7 - 4J=e + 6.5 


-0.5 


c< 

^- x y 


in> 


IT) 


ITI) 


5e + 2Ae + 10(7 - 20J 


e + Ae + 4(7 - 7J=e + 7.5 


0.5 


'-yz 


IU> 


1 Tl) 


IT) 


5e + Ae + 10(7 - 20J 


e + 4(7 - 7J=e + 7.4 


0.4 



corresponds to c yz \GS), and the broad structure at E w —1.3 eV to a superposition of c]. z \GS) and clJGS 1 ). The 
peak at E « eV finally can be identified with ct |G_>). 

Note that the peak in the xy-DOS at E as —0.8 eV has no correspondence in the atomic scheme. It could, however, 
be due to a finite admixture of the state | T)| T)l Tl) to the ground-state of the interacting system. This state would 
then allow for an excitation involving c\, y with an excitation energy —0.6 eV. 
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